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ABSTRACT 

Forces Due to Waves on Submerged Structures. (May 1970) 
George Edward Shank, B.S.C.E., Duke University 
Directed by: Dr. John B. Herbich 

Forces due to gravity waves on several models of half- 
cylindrical and rectangular shape were determined experimentally 
in a two dimensional wave tank. The parameters of water depth, 
wave length and wave height were varied and their effect on forces 
determined. The horizontal and vertical forces on the model were 
measured by instrumented load cells. They were recorded on 
electronic recorders. The results were analyzed with the aid of a 
digital computer. Forces computed using existing theories were 
compared to measured forces and the results used to determine 
inertial coefficient. The results of the data were plotted in 
dimensionless form to provide a correlation between ratios of wave 
length/water depth, wave height/wave length, and dimensionless 
force. The dimensionless force curves can be used to estimate 
forces due to waves on a prototype. 
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INTRODUCTION 

In recent years, offshore oil exoloration and production have 
moved to remote areas of the world far from natural or artificial 
harbors. Sometimes the harbors in the remote areas do not have 
adequate handling capacities to accommodate even small tankers. 

It appears that savings in both capital investment and operating 
costs may be realized if the oil could be stored in submerged 
tanks in the oil producing areas. This conclusion is reached 
by Chamberlin (2) in his consideration of underwater tank designs. 
Also in cases of small, undeveloped offshore fields it may not 
be economical to pump oil using a pipeline to shore and to develop 
shore facilities. 

One of such underwater storage tanks was recently installed in 
the Middle East by the Chicago Bridge and Iron Company. This 
tank has a 500,000 barrel capacity, is about 260 feet in diameter 
and about 65 feet high. It was installed in water 156-feet deep. 
Since such a tank involves a large investment of money and since 
no adequate theoretical solutions are available to date, extensive 
model wave tank investigations were conducted prior to design 
and construction of the tank by the Chicago Bridge and Iron 
Company. 

The Institute of Ocean Development (8) in Tokyo conducted field 

The citations on the following pages follow the style of the 
Journal of the Hydraulics Division , Proceedings of the American 
Society of Civil Engineers. 






studies in Chiba Port adjacent to Tokyo with a model one-tenth 
the size of the prototype tank. The uooer half of the model was 
hemispherical in shape and utilized flexible material. The lower 
half of the model was elipsoidal in shape and was made from steel. 
The model had a capacity of about 10,000 aallons. 

The United States Commission on Marine Science, Engineering 
and Resources (3, 4) and the Oceanographer of the Navy (14) have 
predicted that underwater facilities will be built at an increasing 
rate in the near future. These will be built either in connection 
with oil and gas production, or in connection with mineral recovery 
from the ocean bottom, or as underwater habitats. Submerged mili¬ 
tary support bases for submersibles are predicted. These support 
bases will eliminate sea state as a cause for the prevention or 
termination of operations to support submersibles. Habitats can 
support scientific explorations and the exploitation of the natural 
resources of the ocean floor. A great need exists to be able 
to predict forces caused by waves on submerged structures. Until 
better theoretical force predictions are available, design criteria 
will have to rely upon model tests. 

On the basis of model tests, this thesis provides design in¬ 
formation about wave forces on a selected number of shapes of 


underwater structures. 







LITERATURE SURVEY 


Forces on a body oscillating in a viscous fluid or on a fixed 
body in an oscillating viscous fluid are due to two causes. These 
are inertia and drag. The drag term is a function of velocity, 
the inertial term a function of acceleration. This observation 
was first made by G. G. Stokes (9) in 1851. 

Forces on piles due to gravity waves were first studied by 
Morison, O'Brien, Johnson and Schaaf (1) in 1950. They set forth 
the so-called, and well known "Morison Equation" which relates 
forces to wave parameters. The force was found to consist of 
two components. These were a drag force proportional to the square 
of the velocity which may be represented by a drag coefficient 
having substantially the same value as for steady flow and a 
virtual mass force proportional to the horizontal component of the 
acceleration force exerted on the mass of displaced water. The 
Morison Equation is: 

F = Fp + F^ = l/2Cp p A|u|u + p V 1^.(1) 

where:^ 

F = Total force 
Fp = Drag force 

Vhe literature survey will use the symbols as they appear in 
the cited literature, and will be redefined for each author. Sym¬ 
bols for other sections of the thesis will be consistent and will 
be defined only the first time they appear. A list of symbols 
is given in Appendix I. 
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F., = Inertial force 
M 

Cp = Drag coefficient 

C., = Inertial coefficient 
M 

A = Cross sectional area of object in direction perpen¬ 
dicular to the wave travel 

V = Volume of object 

0 = Mass density 

u = Horizontal velocity of water particle 

^ = Horizontal acceleration of water particle 

a L 

To verify the theory experimentally, a pile was attached by its 
base to the bottom of a wave channel and restrained by a load cell 
at the top. The laboratory results showed close agreement between 
computed and measured force with the exceotion of some experimental 
scatter. Further work on waves on piles was published by Morison 
(11) in 1951 and Morison, Johnson and O'Brien (12) in 1954 
utilizing the same basic relationship. 

In 1952 O'Brien and Morison (13) applied the Morison Equation 
to a submerged spherical object. The components of force were 
due to both drag and inertia. Forces on the sphere were evaluated 
when it was suspended from a rod above the bottom and when It was 
placed on the bottom by use of a tee made of razor blades. The 
tee was made by attaching tv/o razor blades to the tank bottom. 

The blades were perpendicular to the direction of wave travel and 
parallel to each other. The suspended sphere was used to measure 
horizontal force. The sphere on the tee was used to measure 
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vertical force. In the latter case, the period was held constant 
and the wave height raised in small increments until the sphere 
tipped off the tee. The uplift force was equal to the weight of 
the sphere when movement occured. Drag force and inertial force 
were computed by utilizing the measured force in the cases where 
the crest, trough and nodes of the wave were over the sphere. 

In 1953 Reid and Bretschneider (15) utilized the Morison 
Equation for the computation of wave forces on piles and derived 
formulas for computation of wave forces on large submerged objects. 
The forces on large objects were found to be almost entirely 
inertial. This permitted force due to drag to be neglected in 
force computations. Due to the fact that inertial force is a 
function of particle acceleration, once drag is neglected the 
maximum horizontal force occurs when the wave node is over the 
model center. The inertial force was computed directly from the 
horizontal distribution of pressure. The equations developed by 
Reid are simple, straight-forward and accurate for large objects 
where drag forces are negligible and inertial forces predominant. 
For forces on a large submerged object utilizing small amplitude 
(Airy) wave theory: 

F|^ = Cj^L^Lj ^ [cose.j - cos (e^ + -J—!-)].(2) 

where: 

F|^ = Total horizontal force 
= Inertial coefficient 
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= Length of structure (in direction of wave travel) 

L 2 = Width of structure 

= Height of structure 

K = Pressure factor = . 

cosh Zvu 

L 

d = Water depth 
L = Wave length 
Y = Specific weight 
H = Wave height 

e-j = Phase angle (position) of end of model facing waves 
measured from crest in direction of wave travel 
The maximum value of F for this wave will occur when the 
center of the object is at the position ^ ^ Thus: 


ttL 


1 


fMax = W3 (-r> 


(3) 


A re-analysis of all existing data on wave forces on piles 
was made by Crooke (5) in 1955. His paper reconciled differences 
in previous studies by applying a consistent approach which 
resulted in reasonable results. He applied the Iverson approach 
to wave forces. The analysis was made utilizing only one co¬ 
efficient in the wave force equation. This coefficient was 
dependent both upon velocity and acceleration. 


C = 


fli VL 
/ AL _p 

\ 2 ' y 




(4) 






















-* • 
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where: 


AL 

= 

Iversen's modulus 

V 

VL 



_P 

u 

= 

Reynolds number 


V 

gL 

= 

Froude number 

A 

= 

Particle acceleration 

V 

= 

Particle velocity 

L 

= 

Characteristic length 

P 

= 

Density 

M 

= 

Vi scosity 

g 

= 

Gravity 


Crooke's approach worked well for deep and shallow water waves 
as long as the appropriate wave theory was used. 

In 1958 Keulegan and Carpenter (9) investigated the inertial 
and drag coefficients of cylinders and plates when exposed to 
gravity waves. They found that a remainder function AR could 
be introduced to provide a truer representation of force when con¬ 
sidering the coefficients and Cp as being constant through 
the wave cycle. AR was required because the point values of 
Cjyj and Cp differed from their average values. They also found 

that average values of Cp, and C., could be correlated with point 

^m''' 

values in the Morison Equation by introducing a parameter . 

U = Maximum horizontal particle velocity 
m 


where: 
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T = Wave period 

D = Diameter of cylinder or width of plate 
The Mori son Equation thus modified becomes: 




= {- 


<U 


Do 


sine) 


m 


|cose|cose. (5) 


where: 


p = Density 
0 = 2ti/T 
0 = Phase angle 

When aR is introduced into this equation as a remainder it 
becomes: 


F _ /IT ,- \ /Do - ^ , 

.^2 p ■ M ^m^ ^0“ ■ T |cose|cose + aR . . . . (6) 

P u mU m 

AR is obtained by subtracting the comouted values for 

(f sine) --y |cose|cose 

m 

from observed values for 

It was found that inertial coefficient and drag coefficient varied 

with the diameter of the cylinder. For a given cylinder it 

varied with the values of maximum current U . A strong correla- 

U T 

tion between the period parameterand inertial and drag coef- 
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ficients was found to exist. 

It was found that for the plates the higher values of the drag 
coefficient are obtained at the smaller values of 



D 


and the higher values of inertial coefficient are associated with 
larger values of 

U T 
m 

D 

Brater, McNown and Stair (1) investigated exoerimentally the 
effect of gravity waves on models of submerged barge-like structures 
in 1958. They evaluated a flat plate, a right parallelepiped and 
two parallelepipeds with a somewhat modified shape. The models 
were suspended from rods instrumented with load cells. The 
Mori son Equation was utilized for analysis. 


F= 1/2C;,Ab |u|u.C„p V|i, . . . 
where: 

Cp = Drag coefficient 
A = Cross sectional area of model 
p = Density 

u = Horizontal particle velocity 

= Inertial coefficient 
m 

V = Volume of model 


— = Horizontal 


acceleration of oarticle 


(7) 
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X 

The maximum force on the model structures occurred between ^ 
and ^ = 1/4 


where: 


L = Wave length 

X = Distance from wave crest to location of maximum 
force measured in direction of wave travel 

Drag forces were found to be negligible on the barge. Analysis 
of the results were made utilizing the inertial portion only. The 
analysis of inertial force used pressure to evaluate horizontal 
force. 



( 8 ) 


where: 


P, = Pressure on front of model 


P 2 = Pressure on rear of model 


This expression is substituted in the inertial portion of the 
Mori son Equation to yield 



(9) 


where: 


F. = Portion of force due to inertia 


The equation of motion in the vertical direction is 



( 10 ) 











where: 


n 


= Increment of pressure 

o Z 

w = Specific weight of water 

When Eg. 10 is substituted into the Airy expressions for particle 
acceleration and integrated, the result is 


P = w(-z + y) + 


ttLH cosh 2TT{d+z)/L - cosh 27i(d+y)/L /n\ 

-sTnK 2Vd/E- 


where ; 

y = Elevation of surface above still water 
z = Distance from water surface to model 
H = Wave height 
T = Wave period 
d = Depth of water 

For the case of maximum force which occurs when the model is 
centered at the quarter point of the wave, the above equation be¬ 
comes : 

2ttXi p ttH cos 2ttXi/L 

- P 2 = wH cos -j—!- [K - 2 sinh^ {- 

where: 

x^ = Horizontal distance from wave crest to leading 

edge of barge 

y _ cosh 2T((d+z)/L 
cosh 27rd/L 

When measured values are substituted in Eq. 12 the value of 

(P-l - P 2 ) can be computed. Eq. 9 can be rewritten as 
F. 

r - 1 

^m " (P| - p2)Ax 


( 13 ) 
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from which C can be computed, 
m 

Both drag force and inertial force were found to be significant 
in the results from the flat plate. Values of Cp and which 
yielded the best correspondence with measured forces and phase 
angles were utilized in the Morison Equation. The values for 
computed wave parameters were derived utilizing Airy theory. For 
any given wave, the force is a function of the phase angle only. 

The Morison Equation was modified as follows to derive the phase 
angle for maximum force where drag and inertial forces were both 
significant, 

F = Cr,' Cos^ e + C ' Sin 6.(14) 

U m 

where: 

F = Total force 
e = Phase angle 


C 


D 


Cp A p 


2 2 
2 

2T 


r cosh 2tt (d+z)/L- |2 
*- sinh 2ir d/L-^ 


(15) 


m 


= p V 


2tt^H 


m 


[- cosh 2tt (d+z)/Li 

*■ sinh 2it d/L ^ 


(16) 


The phase angle for maximum force is obtained by setting dF/de = 0, 
with the result 


Sin 0 


m 


max 


2C 


D 


( 17 ) 


They found that inertial coefficient for the barge was 1.5. 
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Forces obtained from tests of a barge with rounded edges and a 

barge with a slot in the center were not substantially different 

from those obtained from tests on a normal barge with square 

edges. This was due to the fact that inertial force is a function 

of volume. Measured inertial forces on the barge agrees well with 

theory. The drag coefficient for the plate was 3.5 and virtual 

mass (the product of volume and C ) was 1.75 times that of a 

m 

circumscribing cylinder. 

A recent paper by Evans (6) published results of study of 
forces on piles of full-size prototype offshore platforms. He 
utilized the Mori son Equation to compare computed predictions with 
measured results. The measured data included waves of large size 
(up to 106 ft.) due to hurricanes. His data was taken during two 
extended periods. They were 1954 to 1958 and 1960 to 1963. The 
study proved the accuracy of the Mori son Equation, The average 
values of measured and theoretical forces were within 15 percent 
of each other. The value obtained for was 0.5. The value 
for was 1.5. This agrees well with previous work done by 
others and summarized by Wilson and Reid (17). The drag coef¬ 
ficients were found to decrease for increasino wave height. The 
inertial coefficient did not change greatly with moderate changes 
in wave size. However, it doubled in value from a 23 ft. wave to 


a 106 ft, wave. 



* m 


i?4 .C.- 


■-«* 


■* 4 ^^ ^ 


■i- 
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THFnRFTICAL CnNSinFRATinws 

Progressive linear wave (Airy) theory is used for this study. 

A brief presentation is made of the portions aoprooriate to this 
work. Detailed development is made by Faaleson and Dean (7) and 
Wiegel (16). 

Airy theory makes the following assumptions: 

That the bottom is smooth and impermeable 

That the fluid is incompressible, irrotational and inviscid 

That the depth is constant (measured from water level) 

That surface tension is neglected 
That Coriolis force is neglected 
That atmospheric pressure is constant 

That the height of the surface is small in comparison with 
wave lenath 

That all non-linear terms may he neglected 

The wave surface is described by the following relationship: 

n = ^ cos 2tt (^ - y). (Ifi) 

where: 

n = Distance of the free surface from mean water level 
X = Distance from the wave crest measured in the direction 
of wave travel 
L = Wave length 
H = Wave height 


t = Time 
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T = Wave period 

Some of the various wave measurements and parameters which will 
be used in this section and the results and conclusions section 
are shown in Fig. 1. 

The velocity potential for the small amplitude wave is given 
by the expression 

♦ = ^ cos (kx - «t).(,9) 


where; 

<}' = Velocity potential 
g = Gravity 

a = Wave amplitude = H/2 
a = 2 tt/T 

z = Distance from still water level 
d = Water depth 
k = 2 it/L 

The celerity of the wave is given by the expression: 


C = f. tanh ^ 
dTT L 


[5^ tanh 


2TTd-|l/2 

L 


( 20 ) 


where: 


C = Celerity 

The expression for wave length is; 

L = -^ tanh ^.(21) 

2tt L 

In deep water the depth becomes large relative to wave length 

o J 

and the hyperbolic tangent function (tanh -^) approaches unity. 











I 
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WAVE DEFINITION SKETCH 
FOR AIRY WAVE THEORY 
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Wave celerity becomes 



( 22 ) 


where: 


Cq = Deep water celerity and wave length becomes 



(23) 


where: 


Lq = Deep water wave length 


The wave period is an independent variable. It also remains 
relatively constant irregardless of changes in wave height, water 
depth, and wave length. If corresponding expressions for general 
waves and deep water waves are solved using a constant depth, the 
following expression results 



(24) 


0 


This enables direct solution for the wave length from the deep 
water wave length which is easily calculated from Eg. 23. 

The horizontal component of water particle velocity is given 
by the expression. 



(25) 


where: 


u = Horizontal velocity of water particle 


The horizontal component of water particle acceleration is 











4 


Ml 
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given by the expression; 


u 


3u _ cosh k(z+d) 
Jt - 9^'^ cosh kd 


sin (kx - at) 


(26) 


where: 

u = Horizontal acceleration of water particle 
The Morison Equation relates force to particle motion. 


= 1/2 p Cp A|u|u + p V u. (27) 

where: 

p = Density 

= Horizontal force 

For a given vertical increment of the model, this expression be¬ 
comes 


dFh . dFf, + dF, 


(28) 


where: 


Fp = Force due to drag 
Fj = Force due to inertia 


dF. = C p W L u dz 
I m mm 


(29) 


where: 


W = Width of model 
m 

= Length of model 

H^ = Weight of model 
m ^ 

The total inertial force can be obtained from substitution of 
Eq. 26, in Eq. 29 and integration of the resultant expression. 
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-d + H 


m 


F. = C p W L 
I m ^ m m 


ZiT^h cosh l<(d+ 2 ) 4 . /->n\ 

- -T si'rh kd -5"’ 


2 

(C pW L ) (22^-ii-i2IL£i) 
^ ^ ^ sinh kd 


»-d + H 


m 


cosh k(d+z) dz 


ir II I sin at\/Sinh k(d+z) 

(CmpWmLm^ ^~~2 -^ ^-iT^- 

r sinh kd ^ 


-d + H 


m 


0 2 M -4„ 4. sinh k(d-d+H ) - sinh k(d-d) 

(C pW L )( -?; ” . ?t )(-!!!_- 

" """ f sinh kd ^ 


o i 4 . sinh dk H 
f II j H Sin ot\/ m\ 

^ ^~2 -^ (-k-^ 

r sinh kd 


(31) 


The force due to drag on a vertical increment of model is 


dFf 3 = 1/2 Cp p Wj^luju dz 


(32) 


The total drag force can be obtained by 
Eq. 32 and integration of the resultant 



substitution of Eq. 25 in 
expression. 

|u|u dz.(33) 
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= 1/2 C. p W, 


m 



-d + H 


m 


,cosh k(d+z) j. 


1/2 C. p (lcpi_2licos^^ 

^ "" ' (sinh kd)'^ 


-d + H 


m 


cosh*^ k(d+z) dz 


-d 


1/2 C p W (—)^ ^ |cos at|cos gt ^j- sinh 2k(d+z) ^ k(d+z)j 
^ (sinh kd)^ ^ ^ 


-d + H 


m 


MZ Cr, p W (Jc os ot | c ^ a _ t )( 

(s1nh Mr 


sinh 2k(d-d+H ) k(d-d+H ) 


in 


m 


sinh 2k(d-d) k(d-d) n 
4 ~ 2 ’ 


1/2 Cp, p W ( Ic o s . g . t | cos 

^ ^ (sinh kd)^ 


sinh 2k H„ k H 


-) . . . (34) 


If a submerged object has a large volume, the inertial forces 
predominate, and the drag forces can be ignored. If it has a 
length in the direction of wave travel which is large compared to 
the wave length, an analysis of force can be computed directly 
from the horizontal distribution of pressure beneath the wave. 






















The difference between the actual gauge pressure at the bottom 
and the calculated hydrostatic pressure based on the still water 
depth is defined as pressure anomaly. 

aP = P - yd.(35) 

where: 

aP = Pressure anomaly 

P = Gauge pressure 

Y = Specific weight of water 

A pressure factor which is a function of H/d is defined as follows 


(ap)^ - (ap)^ 

7h 


(36) 


where: 

K = Pressure factor 

(ap)^ = Pressure anomaly under the wave crest 
(ap)^ = Pressure anomaly under the wave trough 


. /27rd\ 

cosh (-j^) 


(37) 


For Airy wave theory the variation of Ap/yH with phase angle is a 

|/ 

function (^) cos 6 where e is the phase angle, and is equal to 
at or kx. 

The general equation for horizontal force acting on a large 


rectangular submerged object is 









t 
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where: 

aP^ = Pressure anomaly at leading edge of object 
AP 2 = Pressure anomaly at trailing edge of object 

In the special case of the Airy theory wave 


Fh = ^ [cose, - cos(e, + 2tt 

h mmm2 I 1 L 


L 




(39) 


where: 


0 ^ = ^""^l 


L 

= Horizontal distance from the leadina edge in direc¬ 
tion of wave travel relative to a crest 

The maximum value of the horizontal force (Fn) , occurs when the 

H max 

center of the tank is at the position 

ft = ^ 

0 2 

radians. Therefore, 




( 40 ) 
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EXPERIMENTAL EQUIPMENT AND PROCEDURE 


Models 

With the exception of the flat plate the models were all 
constructed from plexiglass. The flat plate was fabricated 
from aluminum. All models used (with the exception of the aluminum 
plate) are shown in Fig. 2 and Fig. 3. The models consisted of 

(a) plate 7-23/32 in. wide, 4-1/2 in. high and 5/8 in. thick. 

(b) rectangular tank 7-23/32 in. wide, 4-1/2 in. high and 2 in. 
long. 

(c) rectangular tank 7-23/32 in. wide, 4-1/2 in. high and 4 
in. long. 

(d) rectangular tank 7-23/32 in. wide, 4-1/2 in. high and 
8 in. long. 

(e) a half-cylinder tank having radius equal to 4-15/32 in. 
and length equal to 7-11/16 in. 

Equipment and Procedure 

The models were placed in a three feet deep, two feet wide and 
one hundred and twenty feet long wave channel. Waves were generated 
with an oscillating pendulum-type generator capable of producing 
a wide range of wave lengths and heights. The generator is powered 
with a variable speed motor. The wave period of the waves was 
varied by changing the speed of the motor. The wave generator 
is shown in Fig. 4. 

The wave height was changed by varying the adjustment on the 
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FIG. 2.-MnnFLS STuniEn 




FIG. 3.-CLOSE-UP OF FOUR INCH LONG RECTANGULAR MOOFL 
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FIG. 4.-WAVE GENERATOR 



FIG. 5.-EXPERIMENTAL FACILITY 
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excentricity of the paddle arm on the flywheel which resulted in 
a change in the stroke of the pendulum. The pendulum was adjusted 
to provide the proper shape for each wave. A stroke which had 
equal amplitude from top to bottom was used for the most shallow 
water waves. A stroke with a large amplitude at the top and a small 
amplitude at the bottom was used for deep water waves. An efficient 
absorber is installed at the other extremity of the tank to prevent 
reflected waves from interfering v/ith generated incident waves. 

The absorber has a permeable face with an impermeable backina. The 
channel has 3/8 inch glass side panels throughout its length to 
allow for visual observation of experiments. 

The experiments were conducted at v/ater depths of 24 inches, 

18 inches and 13 inches. Wave lengths were varied to produce 
several L/d ratios between 1.5 and 6.0. This covers deep water 
waves and the deep half of intermediate waves (where intermediate 
waves are defined as having a L/d ratio between 2 and 20). The 
pendulum was timed by stopwatch through 20 to 50 cycles and the 
period determined to the closest 1/100 second. 

Approximately 100 different combinations of water depth, wave 
period (and consequently wave length) and wave height were used for 
each model. The wave heights were increased until the waves broke. 
The schedule of test waves used are shown in Table 1. For each 
combination of water depth and wave period the schedule increases 
wave height to a value which will yield a steepness equal to the 
theoretical maximum of 1.43, A total of 599 different experimental 
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TABLE 1. - TEST WAVE SCHEDULE 


Test 

Wave 

Number 

Depth 
i n 

feet 

Period 
in 

seconds 

Wave 

Hei ght 
in feet 

d/L 

L/d 

H/L 

1 

2.000 

1.68 

0.10 

0.178 

5.61 

0.009 

2 

2.000 

1.68 

0.20 

0.178 

5.61 

0.018 

3 

2.000 

1.68 

0.30 

0.178 

5.61 

0.027 

4 

2.000 

1.68 

0.40 

0.178 

5.61 

0.036 

5 

2.000 

1.68 

0.50 

0.178 

5.61 

0.044 

6 

2.000 

1.68 

0.60 

0.178 

5.61 

0.053 

7 

2.000 

1.68 

0.70 

0.178 

5.61 

0.062 

8 

2.000 

1.55 

0.10 

0.203 

4.93 

0.010 

9 

2.000 

1 .55 

0.20 

0.203 

4.93 

0.020 

10 

2.000 

1.55 

0.30 

0.203 

4.93 

0.030 

11 

2.000 

1.55 

0.40 

0.203 

4.93 

0.041 

12 

2.000 

1.55 

0.50 

0.203 

4.93 

0.051 

13 

2.000 

1.55 

0.60 

0.203 

4.93 

0.061 

14 

2.000 

1.55 

0.70 

0.203 

4.93 

0.071 

15 

2.000 

1 .29 

0.10 

0.265 

3.77 

0.013 

16 

2.000 

1.29 

0.20 

0.265 

3.77 

0.027 

17 

2.000 

1.29 

0.30 

0.265 

3.77 

0.040 

18 

2.000 

1.29 

0.40 

0.265 

3.77 

0.053 

19 

2.000 

1.29 

0.50 

0.265 

3.77 

0.066 

20 

2.000 

1.29 

0.60 

0.265 

3.77 

0.080 

21 

2.000 

1.29 

0.70 

0.265 

3.77 

0.092 

22 

2.000 

1.09 

0.10 

0.342 

2.92 

0.017 

23 

2.000 

1.09 

0.20 

0.342 

2.92 

0.034 

24 

2.000 

1.09 

0.30 

0.342 

2.92 

0.051 

25 

2.000 

1 .09 

0.40 

0.342 

2.92 

0.068 

26 

2.000 

1.09 

0.50 

0.342 

2.92 

0.086 

27 

2.000 

1.09 

0.60 

0.342 

2.92 

0.103 

28 

2.000 

0.89 

0.10 

0.505 

1.98 

0.025 

29 

2.000 

0.89 

0.20 

0.505 

1 .98 

0.051 

30 

2.000 

0.89 

0.30 

0.505 

1.98 

0.076 

31 

2.000 

0.89 

0.40 

0.505 

1.98 

0.101 

32 

2.000 

0.89 

0.50 

0.505 

1.98 

0.126 
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TABLE 1. - TEST WAVE SCHEDULE 
(Continued) 


Test 

Wave 

Number 

Depth 

in 

feet 

Period 
i n 

seconds 

Wave 
Height 
in feet 

d/L 

L/d 

H/L 

33 

2.000 

0.78 

0.10 

0.646 

1 .55 

0.032 

34 

2.000 

0.78 

0.20 

0.646 

1.55 

0.065 

35 

2.000 

0.78 

0.30 

0.646 

1.55 

0.097 

36 

2.000 

0.78 

0.40 

0.646 

1.55 

0.129 

37 

2.000 

0.78 

0.50 

0.646 

1.55 

0.161 

38 

1.500 

0.71 

0.10 

0.579 

1.73 

0.039 

39 

1.500 

0.71 

0.20 

0.579 

1.73 

0.077 

40 

1.500 

0.71 

0.30 

0.579 

1.73 

0.116 

41 

1.500 

0.71 

0.40 

0.579 

1.73 

0.154 

42 

1.500 

0.79 

0.10 

0.467 

2.14 

0.031 

43 

1.500 

0.79 

0.20 

0.467 

2.14 

0.062 

44 

1.500 

0.79 

0.30 

0.467 

2.14 

0.093 

45 

1.500 

0.79 

0.40 

0.467 

2.14 

0.125 

46 

1.500 

0.79 

0.50 

0.467 

2.14 

0.156 

47 

1.500 

0.92 

0.10 

0.353 

2.84 

0.023 

48 

1.500 

0.92 

0.20 

0.353 

2.84 

0.047 

49 

1.500 

0.92 

0.30 

0.353 

2.84 

0.070 

50 

1.500 

0.92 

0.40 

0.353 

2.84 

0.094 

51 

1.500 

0.92 

0.50 

0.353 

2.84 

0.117 

52 

1.500 

0.92 

0.60 

0.353 

2.84 

0.140 

53 

1 .500 

1.09 

0.10 

0.265 

3.78 

0.018 

54 

1.500 

1.09 

0.20 

0.265 

3.78 

0.035 

55 

1.500 

1.09 

0.30 

0.265 

3.78 

0.053 

56 

1.500 

1.09 

0.40 

0.265 

3.78 

0.071 

57 

1.500 

1.09 

0.50 

0.265 

3.78 

0.088 

58 

1.500 

1 .09 

0.60 

0.265 

3.78 

0.106 

59 

1.500 

1.31 

0.10 

0.200 

5.01 

0.013 

60 

1.500 

1 .31 

0.20 

0.200 

5.01 

0.027 

61 

1.500 

1.31 

0.30 

0.200 

5.01 

0.040 

62 

1.500 

1 .31 

0.40 

0.200 

5.01 

0.053 

63 

1 .500 

1.31 

0.50 

0.200 

5.01 

0.067 

64 

1.500 

1.31 

0.60 

0.200 

5.01 

0.080 

65 

1.500 

1.31 

0.70 

0.200 

5.01 

0.093 
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TABLE 1. - TEST WAVE SCHEDULE 
(Continued) 


Test 

Wave 

Number 

Depth 
i n 

feet 

Peri od 
in 

seconds 

Wave 

Height 
in feet 

d/L 

L/d 

H/L 

66 

1.500 

1.47 

0.10 

0.171 

5.85 

0.011 

67 

1.500 

1.47 

0.20 

0.171 

5.85 

0.023 

68 

1.500 

1.47 

0.30 

0.171 

5.85 

0.034 

69 

1.500 

1.47 

0.40 

0.171 

5.85 

0.046 

70 

1.500 

1.47 

0.50 

0.171 

5.85 

0.060 

71 

1.500 

1.47 

0.60 

0.171 

5.85 

0.068 

72 

1.500 

1.47 

0.70 

0.171 

5.85 

0.080 

73 

1.083 

1.26 

0.10 

0.169 

5.93 

0.016 

74 

1.083 

1.26 

0.20 

0.169 

5.93 

0.031 

75 

1.083 

1.26 

0.30 

0.169 

5.93 

0.047 

76 

1 .083 

1.26 

0.40 

0.169 

5.93 

0.062 

77 

1.083 

1.26 

0.50 

0.169 

5.93 

0.078 

78 

1 .083 

1.26 

0.60 

0.169 

5.93 

0.093 

79 

1.083 

1.04 

0.10 

0.220 

4.54 

0.620 

80 

1.083 

1.04 

0.20 

0.220 

4.54 

0.041 

81 

1.083 

1.04 

0.30 

0.220 

4.54 

0.061 

82 

1.083 

1.04 

0.40 

0.220 

4.54 

0.081 

83 

1 .083 

1.04 

0.50 

0.220 

4.54 

0.102 

84 

1.083 

1.04 

0.60 

0.220 

4.54 

0.122 

85 

1 .083 

0.89 

0.10 

0.283 

3.53 

0.026 

86 

1.083 

0.89 

0.20 

0.283 

3.53 

0.052 

87 

1 .083 

0.89 

0.30 

0.283 

3.53 

0.078 

88 

1 .083 

0.89 

0.40 

0.283 

3.53 

0.105 

89 

1 .083 

0.89 

0.50 

0.283 

3.53 

0.131 

90 

1.083 

0.89 

0.60 

0.283 

3.53 

0.157 

91 

1.083 

0.75 

0.10 

0.382 

2.62 

0.035 

92 

1.083 

0.75 

0.20 

0.382 

2.62 

0.070 

93 

1 .083 

0.75 

0.30 

0.382 

2.62 

0.106 

94 

1.083 

0.75 

0.40 

0.382 

2.62 

0.141 

95 

1.083 

0.71 

0.10 

0.421 

2.37 

0.039 

96 

1.083 

0.71 

0.20 

0.421 

2.37 

0.078 

97 

1.083 

0.71 

0.30 

0.421 

2.37 

0.117 

98 

1.083 

0.71 

0.40 

0.421 

2.37 

0.156 
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runs were made. The resulting massive amount of data provided 
almost total coverage within the range of wave parameters tested. 

An over-all view of the experimental facility is shown in 
Fig. 5. A schematic diagram of the experimental facility is shown 
in Fig. 6. To accomplish the force measurements the model was 
supported by four vertical wires with the upper ends of the wires 
supported from cantilevered load cells. The wire was 0.015 inch 
diameter stainless steel cable. It was composed of 7 strands and 
provided a supple, inelastic support for the model. It was 
fastened to the model by brass and stainless steel swivel snap 
hooks at the lower end and connected to the load cell by a threaded 
eye bolt at the upper end. The cable was fastened to the snap 
hook and threaded eye by use of a crimped stainless steel sleeve. 
The threaded eye had 1 1/2 inches of threaded shaft to facilitate 
vertical adjustment of the models. The models were placed with 
their base 1/4 inch above the tank bottom. A close-up view of 
the experimental set up is shown in Fig. 7. 

Vertical forces caused the load cells to move in a vertical 
direction. The load was measured by strain gauges and recorded 
on electronic recorders. 

Since wave forces produce loads both ud and down on the model, 
it was necessary to have it weighted with metal inside the model 
so that the vertical load lines were pretensioned. This pre¬ 
tension was of high enough value so that the lines never went 
slack. By hanging the structure from four load lines it was also 
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FIG. 7.-CLOSE-UP OF EXPERIMENTAL SET-UP 
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possible to measure moments on the model. A model suspended from 
the vertical load cells is shown in Fig. 8. 

The horizontal forces were measured as follows: A wire and 
spring were used to attach the forward end of the model to the wave 
tank bottom. The downstream end of the model was attached by a 
wire which was diverted ninety degrees over a pulley to a load cell. 
The spring provided a constant tension, the variation of which was 
measured by the load cell and indicated the wave induced horizontal 
force on the model. The measured load was recorded on an electronic 
recorder. 

Both the stud and the pulley were placed at a horizontal dis¬ 
tance of three feet from the model to minimize disturbance in the 
model area. A large diameter (5 inches) ball bearing pulley was 
used to minimize the effects of friction. A schematic diagram of 
the model and load cells is shown in Fig. 9. Known weights were 
used to calibrate the force measurement apparatus. A capacitance 
type wave gauge was used to measure wave profiles. The wave gauge 
was calibrated by immersing the gauge wire into the water by incre¬ 
mented known distances with the water still in the channel. The 
wave height and wave profiles were recorded on an electronic 
recorder. 

A probe was positioned adjacent to the model so that the 
passing of each wave crest past the model would be marked on the 
recorder. This made possible the measurement of the phase angle 
(or time difference) between the crest and the maximum wave forces. 
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FIG. 8.-HALF-CYLINDER MODEL IN PLACE 
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A synchronized marker was used to place a common time mark simul¬ 
taneously on each recorder data paper. This enabled the data 
outputs from the three recorders to be correlated in time. The 
data was recorded on hot needle-type direct-recording equipment. 

The wave height and wave profiles were recorded on a Hewlett-Packard 
dual-channel carrier pre-amplifier recorder (model 321). The 
wave forces were recorded on Sanborn two-and four-channel recorders 
(model 150). Recorder outputs are shown in Fig. 10 and Fig. 11. 


V' 
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FIG. 10.-RECORDER OUTPUT FOR VERTICAL FORCES AND WAVE CREST INDICATOR 
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FIG. 11.-RECORDER OUTPUT FOR HORIZONTAL FORCES, VERTICAL FORCES AND WAVE PROFILE 
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PRESENTATION AND DISCUSSION OF RESULTS 

The primary objective of the research was to determine ex¬ 
perimentally the forces due to gravity waves on several submerged 
structures of rectangular and half cylindrical shape. The relation 
between measured vertical and horizontal force and various wave 
parameters was investigated. Forces determined by experiment were 
compared to those calculated using existing theories for wave 
forces on submerged objects. Wave reflections from the model were 
assumed to be small. Therefore, they were not measured and were 
neglected in the theoretical analysis. 

The raw data was in the form of recorder output shown in 
Fig. 10 and Fig. 11. The period was determined by measurement 
of the length in decimeters between a wave crest, and the crest 
of the 20th wave following it. The paper feed was 5 millimeters 
per second, so twenty waves allowed reading the measured decimeter 
length as a wave period in seconds. The period thus determined was 
compared to the period measured by timing the pendulum. 

An average of at least 10 waves ms taken for the maximum 
wave height, maximum horizontal force and each of the four 
vertical maximum vertical forces (one for each supporting wire). 

The upward force for each wire was measured separately from its 
downward force. The force in the direction of wave travel was 
measured separately from the force in the opposite direction of 
wave travel. 

Once the raw data was reduced to numerical values, it was 


i 
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reduced by use of an IBM model 360/65 digital computer. The inputs 
to the computer program were run serial number, wave period, water 
depth, wave height, measured horizontal force in the direction of 
wave travel, measured horizontal force in the direction opposite 
from wave travel, measured vertical force down, vertical force up 
and phase angle between the wave crest and maximum forces. The 
program averaged the wave period for each grouo of waves where the 
water depth was constant and the period nominally constant. From 
the average period, the water depth and the wave height, the 
program computed the depth to length (d/L) ratio, the length to 
depth (L/d) ratio and the steepness (H/L). 

The units of the measured forces were grams. These measured 
forces were converted to dimensionless form by the program. The 
dimensionless force was calculated as follows: 


''dim 


F 

^ max 


^ H 
^ d 2 


where: 


(41) 


FpiM = Dimensionless force 

F = Maximum force measured 
max 

A = Significant linear dimension equal to height of the 
model squared divided by length of model in 
direction of wave propagation for all models except 
flat olate 

A = H = Height of model for flat plate 
m ^ 
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The computer program which calculates wave parameters and dimension¬ 
less forces and a sample of computer output are presented in 
Appendix VI. 

A dimensionless plot was made of dimensionless force versus 
steepness for each group of waves with constant depth and period 
(constant L/d ratio). A sample plot of dimensionless force versus 
steepness is shown in Fig. 12. The horizontal force with the 
direction of the wave v/as usually slightly hioher than the hori¬ 
zontal force exerted against the direction of the wave. As the 
forces in the two directions closely approximated each other in 
magnitude, a single curve was used to represent both sets of points. 

A cross-plot was made from the dimensionless force versus 
steepness plot. This cross-plot summarized data by plotting 
dimensionless force versus L/d ratio for wave steepnesses of 0.02 to 
0.07. These summary plots are made for vertical and horizontal 
forces for all models with the exception of the flat plate. Only 
horizontal forces were plotted for the flat plate. These curves 
are presented in Appendix II. Since data are in dimensionless 
form, they can be extrapolated to any size rectangular object or 
to any size of a geometrically similar half-cylinder object for 
any wave height, wave length and water depth. An illustration of 
the use of the curves is given in Appendix V. 

An attempt was made to evaluate both inertial and drag 
coefficient by use of the Mori son Equation. This evaluation was 
made by computer solution of Eg. 28, Eq. 31 and Eq. 34. The 
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TEST NUMBERS 
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FIG. 12- SAMPLE DIMENSIONLESS HORIZONTAL 
FORCE PLOT FOR FLAT PLATE. 
L/d = 5,61 












































43 


inputs to this solution were the same as those cited above for 

dimensionless force computation. The proaram increased the 

values of C. and C by iteration. This was done by iterating 
dm 

larger values of C^, then holding them constant, and comparing 
the total measured force and the computed total force for in¬ 
creasing values of C^. When the measured force and computed force 
approximate each other, the comouted value of is given. In 
the case where drag force is significant the best values of 
and C can be selected by picking the values which remain more 
or less constant with changing wave parameters. For the flat 
plate and two inch long rectangular model, the volume of a cir¬ 
cumscribing cylinder was used for the computations. It was found 
that drag forces were small in the case of the 5/8 inch thick 
flat plate. A value of 3.5 for was suggested by Brater, McNown 
and Stair (1). Using this value, the drag force represented 
less than 5 percent of the total force for waves of small steep¬ 
ness and 10 percent of the total force for waves of large steep¬ 
ness. Apparently, the finite thickness of the plate was enough 
to make inertial force predominant. Due to the small size of the 
drag forces it was impossible to evaluate drag coefficient. For 
all models except the flat plate, the drag force was insignificant. 
The computer program which applies the Mori son Eguation to deter¬ 
mine drag and inertial coefficients and a sample of computer 
output are presented in Appendix VIII. 

Since forces on all the models are largely inertial, the 
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values of coefficient of inertia were computed following the 

methods suggested by Reid and Bretschneider (15). A computer 

program was used which utilized the same inputs as those cited 

above for dimensionless force computation. The program determined 

C by use of Eq. 40. Values of were increased by iteration 
m m 

until measured force and the computed force aoproximated each 

other. This method was used to determine inertial coefficients 

for all models except the flat plate. The computer program which 

computes utilizing Eq. 40 and a sample of computer output 

are presented in Appendix VI. The subprogram titled "Force-2" 

performs this function. The program was modified to sum a series 

of thin rectangular solids to accomodate the half cylinder model, 

as it was not of constant cross section in the direction of wave 

travel. This program can be applied to any large submerged object 

of irregular cross section. The program thus modified and a sample 

of computer output are presented in Appendix VII. The subprogram 

"Force-4" performs this function. The results of calculation of 

by Eq. 40 are presented graphically in Appendix III. These 

curves indicate that the coefficients of inertia show a small 

variation with L/d and H/L ratios. In general the coefficient 

was higher for higher values of H/L. There was also some 

variation of C with L/d as indicated in the curves. The values 
m 

of coefficient C varied between 1.20 and 2.20 for all tests, 
m 

Methods suggested by Reid and Bretschneider (15) were used to 
compare the computed and experimental values. Several values of 
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C were selected for this comparison. The computation utilized 
m 

Eq. 40 to determine the calculated force for a given C . Sub- 

program "Force-1" performed this computation in the computer 

program for objects of regular cross section. This program is 

presented in Appendix VI. Subprogram "Force-3" performed this 

computation in the computer program for objects of irregular cross 

section. This program is presented in Appendix VII. The results 

of the comparison of measured force and calculated force are 

presented in tabular form in Appendix IV. The best correlation 

of the computed and experimental data for the two-inch long model 

was obtained with a value of C of 1.40. The best correlation 

m 

for all other models was obtained with a value of C of 1.80. 

m 


I 
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CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER RESEARCH 

It was found that force was a regular function of steepness 
(H/L) and relative wave length (L/d). The force increased for 
increasing values of either parameter. The general shape was the 
same for all curves of steepness versus dimensionless force. 

The shape of all relative wave length versus dimensionless force 
was similar. As the data is in dimensionless form, it can be 
extrapolated to prototypes of similar shape. 

The data was analyzed utilizina the Mori son Equation. It 
was found that drag force was small for a flat plate of finite 
thickness and negligible for models with larger volumes. The 
inertial portion of the Morison Equation and equations derived 
from it were used for analysis of forces and determination of 
inertial coefficient. The coefficients of inertia varied little 
with changes in relative wave length and steepness. In general 
the coefficient was slightly higher for higher steepnesses. The 
values of inertial coefficient varied between 1.20 and 2.20 
for all tests. Values of C^^ in this range were used to compute 
force. This force was compared to the measured values. The best 
correlation for the shortest (in the direction of wave travel) 
model was obtained with a value of of 1.40. The best cor¬ 
relation for all other models was obtained with a value of of 
1.80. 

As the flat plate tested yielded primarily inertial forces. 


• *^1 


gii 
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future experiments should include tests of a thinner plate. This 
would facilitate evaluation of drag coefficient. 

Prototypes are exposed to the irregular gravity waves of the 
sea, and not the regular ones of the laboratory. For this reason, 
future research should include tests of models exposed to irregular, 
wind-generated waves. These tests could be related to existinq 
work by techniques of spectral analysis. 

It is also suggested that future tests be conducted in a 
three-dimensional wave tank to eliminate any possible errors 
introduced by reflections induced by the walls of a two-dimensional 
wave channel. 

The most valuable future research would be the instrumentation 
of prototypes and very large models. This has already provided 
valuable information in the case of offshore platforms. While 
it is a difficult and expensive proposal, only when prototypes 
are tested will absolute information about forces on submerged 


structures be known. 
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APPRNniX I 

NOTATION 


The following symbols are used in this paper: 


1 


Symbol 

a 

A 


A 

0 

d 

Fr 


DIM 


max 


q 

H 

H 


m 


Quantity 

Wave amplitude = H/2 

Significant linear dimension eoual to 
height of the model squared divided 
by length of model in direction of 
wave prooaqation for all models except 
flat olate 

Height of model for flat plate 

Celerity 

Water depth 

Force due to drag 

Dimensionless force 

Horizontal force 

Force due to inertia 

Maximum force measured 

Gravity 

Wave height 

Height of model 

Pressure coefficient 


Dimension 

L 


L 

L 

L/T 

L 

ML/T 


2 


ML/T 

ML/T 

ML/T 

2 

L/T'^ 

L 

L 


2 

2 

2 


These symbols are consistent thrnuoHout the thesis with the 
exception of the literature survey. In the literature survey, 
symbols are used as they appear in the cited literature, and are 
redefined for each author. 
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Symbol 

k 

L 

Lo 

P 


ap 

ap 


1 


AP 


2 


T 

t 


u 

u 


W 


m 


X 


z 

Y 

n 


Quantity 


2Tr/L 

Wave length 
Lenath of model 
Deep water wave length 
Rauae pressure 
Pressure anomaly 

Pressure anomaly at leading edge of 
object 

Pressure anomaly at trailing edge o-f 
object 

Pressure anomaly under the wave crest 
Pressure anomaly under the wave trough 
Wave oeriod 
T ime 

Horizontal velocity of water oarticle 

Horizontal acceleration of water 
particle 

Width of model 

The distance from the wave crest 
measured in the direction of wave 
travel 

Horizontal distance from the wave crest 
to the edge of the model direction of 
wave travel 

Distance from still water level 

Snecific weight of water 

Distance of the free surface -from 
mean water level 


Dimension 

1/L 

L 

L 

L 

M/LT^ 

M/LT^ 

M/LT^ 

M/LT^ 

M/LT^ 

M/LT^ 

T 

T 

L/T 

L/t2 

L 

L 


L 

L 

m/lV 

L 
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Symbol 


Qiianti tv 


2itXi/L 

Density 

2tt/T 

Velocity potential 


nimension 


M/L' 


1/T 

L^/T 





i' 
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APPFMniX II 

niMENSIONLFSS FORCE SUMMARY CUPVFS 



FIG. 13-DIMENSIONLESS HORIZONTAL FORCE PLOT FOR 
FLAT PLATE. 
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L/d 


FIG. 14-DIMENSIONLESS HORIZONTAL FORCE (COMPUTED AS A 
RECTANGULAR MODEL) FOR 2 INCH RECTANGULAR 
MODEL 
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FIG. 15-DIMENSIONLESS VERTICAL FORCE (COMPUTED AS A 
RECTANGULAR MODEL) FOR 2 INCH RECTANGULAR 
MODEL 
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L/d 

FIG. 16-DIMENSIONLESS HORIZONTAL FORCE (COMPUTED AS 

A FLAT PLATE) FOR 2 INCH RECTANGULAR MODEL 
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>^DIM 



FIG. 17-DIMENSIONLESS VERTICAL FORCE (COMPUTED AS 
A FLAT PLATE) FOR 2 INCH RECTANGULAR 
MODEL 
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L/d 


FIG. 18-DIMENSIONLESS HORIZONTAL FORCE FOR 4 INCH 
RECTANGULAR MODEL 
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'^DIM 



L/d 


FIG. 19-DIMENSIONLESS VERTICAL FORCE FOR 4 INCH 
RECTANGULAR MODEL 
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L/d 

FIG. 20-DIMENSI0NLESS HORIZONTAL FORCE FOR 8 
INCH RECTANGULAR MODEL 
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I^DIM 



L/d 


FIG. 21- DIMENSIONLESS VERTICAL FORCE FOR 
8 INCH RECTANGULAR MODEL 

































L/d 


FIG. 22-DIMENSIONLESS HORIZONTAL FORCE PLOT FOR 
HALF CYLINDER MODEL (END FACING INCIDENT 
WAVE) 
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L/d 

FIG. 23-DIMENSIONLESS VERTICAL FORCE 
PLOT FOR HALF CYLINDER MODEL 
(END FACING INCIDENT WAVE) 



















L/d 


FI6.24-DIMENSI0NLESS HORIZONTAL FORCE PLOT FOR 
HALF CYLINDER MODEL (SIDE FACING INCIDENT 
WAVE) 
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L/d 


FIG.25-DIMENSI0NLESS VERTICAL FORCE PLOT FOR HALF 
CYLINDER MODEL (SIDE FACING INCIDENT WAVE) 
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APPENDIX III 

TNERTIAL CnEFFICIEMT rilPVPS 
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FIG.26-COEFFICIENT OF INERTIA FOR RECTANGU¬ 
LAR MODEL 2-INCHES LONG (WATER 
DEPTH 13 INCHES) 
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FIG. 27-COEFFICIENT OF INERTIA FOR RECTANGU¬ 
LAR MODEL 2-INCHES LONG (WATER 
DEPTH 18 INCHES) 
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FIG. 28-COEFFICIENT OF INERTIA FOR RECTANGU¬ 
LAR MODEL 2 -INCHES LONG (WATER 
24 INCHES) 
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FIG. 29-COEFFICIENT OF INERTIA FOR RECTANGU¬ 
LAR MODEL 4 - INCHES LONG (WATER 
DEPTH 13 INCHES) 
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H/L 


FIG. 30-COEFFICIENT OF INERTIA FOR RECTANGU¬ 
LAR MODEL 4-INCHES LONG (WATER 
DEPTH 18 INCHES) 
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H/L 


FIG. 31-COEFFICIENT OF INERTIA FOR RECTANGU¬ 
LAR MODEL 4-INCHES LONG (WATER 
DEPTH 24 INCHES) 
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FIG. 32- COEFFICIENT OF INERTIA FOR 
RECTANGULAR MODEL 8-INCHES 
LONG (WATER DEPTH 13 INCHES) 
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FIG. 33- COEFFICIENT OF INERTIA FOR 
RECTANGULAR MODEL 8-INCHES 
LONG (WATER DEPTH 18 INCHES) 
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FIG. 34 - COEFFICIENT OF INERTIA FOR 
RECTANGULAR MODEL 8-INCHES 
LONG (WATER DEPTH 24 INCHES) 
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H/L 

FIG. 35-COEFFICIENT OF INERTIA FOR HALF 

CYLINDER MODEL (END FACING INCIDENT 
WAVE, WATER DEPTH 13 INCHES) 
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FIG. 36-COEFFICIENT OF INERTIA FOR HALF 

CYLINDER MODEL (END FACING INCIDENT 
WAVE WATER DEPTH 18 INCHES 
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FIG. 37-COEFFICIENT OF INERTIA FOR HALF 

CYLINDER MODEL (END FACING INCIDENT 
WAVE, WATER DEPTH 24 INCHES) 
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H/L 

FIG. 38-COEFFICIENT OF INERTIA FOR HALF 

CYLINDER MODEL (SIDE FACING INCIDENT 
WAVE, WATER DEPTH 13 INCHES) 
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FIG. 39-COEFICIENT OF INERTIA FOR HALF 

CYLINDER MODEL (SIDE FACING INCIDENT 
WAVE, WATER DEPTH 18 INCHES) 
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FIG.40-COEFFICIENT OF INERTIA FOR HALF 

CYLINDER MODEL (SIDE FACING INCIDENT 
WAVE, WATER DEPTH 24 INCHES) 
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APPENDIX IV 

COMPARISON OF MEASURED AND COMPUTED HORIZONTAL FORCES 


COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR RECTANGULAR MODEL 8 INCHES LONG 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in qrams 

m 

C =1 -8 

m 

C,„=2.0 

m 

2.000 

0.182 

5.496 

0.007 

94 

89 

107 

119 

2.000 

0.182 

5.496 

0.016 

224 

194 

232 

258 

2.000 

0.182 

5.496 

0.023 

366 

276 

331 

368 

2.000 

0.182 

5.496 

0.032 

528 

387 

465 

517 

2.000 

0.182 

5.496 

0.044 

720 

534 

641 

713 

2.000 

0.182 

5.496 

0.056 

850 

690 

828 

920 

2.000 

0.182 

5.496 

0.065 

1000 

791 

949 

1050 

2.000 

0.354 

2.824 

0.030 

100 

132 

159 

176 

2.000 

0.354 

2.824 

0.045 

210 

201 

242 

286 

2.000 

0.354 

2.824 

0.067 

342 

299 

359 

398 

2.000 

0.354 

2.824 

0.092 

460 

409 

491 

545 

2.000 

0.354 

2.824 

0.092 

524 

407 

488 

542 

2.000 

0.354 

2.824 

0.103 

585 

456 

547 

608 

2.000 

0.656 

1.524 

0.026 

20 

17 

20 

22 

2.000 

0.656 

1.524 

0.051 

42 

33 

39 

43 

2.000 

0.656 

1.524 

0.074 

66 

47 

56 

62 

2.000 

0.656 

1 .524 

0.100 

96 

63 

76 

85 

2.000 

0.656 

1.524 

0.112 

106 

71 

85 

95 

2.000 

0.656 

1.524 

0.112 

98 

71 

85 

95 

1.500 

0.577 

1.734 

0.023 

26 

24 

28 

31 

1.500 

0.577 

1 .734 

0.045 

56 

45 

54 

61 

1.500 

0.577 

1 .734 

0.066 

88 

67 

81 

90 

1.500 

0.577 

1 .734 

0.088 

122 

89 

107 

119 

1.500 

0.577 

1.734 

0.111 

146 

113 

136 

151 

1.500 

0.577 

1.734 

0.106 

162 

107 

129 

143 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR RECTANGULAR MODEL 8 INCHES LONG 
(Continued) 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in grams 

C =1.5 
m 

C =1.8 
m 


1.500 

0.171 

5.851 

0.009 

138 

118 

141 

157 

1.500 

0.171 

5.851 

0.021 

310 

267 

321 

357 

1.500 

0.171 

5.851 

0.031 

516 

400 

480 

533 

1.500 

0.171 

5.851 

0.039 

708 

508 

610 

678 

1.500 

0.171 

5.851 

0.055 

920 

711 

853 

948 

1.500 

0.171 

5.851 

0.063 

1080 

811 

973 

1080 

1.500 

0.171 

5.851 

0.074 

1200 

949 

1140 

1270 

1.083 

0.167 

6.002 

0.012 

160 

153 

183 

204 

1.083 

0.167 

6.002 

0.012 

170 

153 

183 

204 

1.083 

0.167 

6.002 

0.012 

192 

155 

186 

206 

1.083 

0.167 

6.002 

0.025 

408 

322 

386 

429 

1.083 

0.167 

6.002 

0.038 

644 

499 

598 

665 

1.083 

0.167 

6.002 

0.051 

952 

668 

801 

890 

1.083 

0.167 

6.002 

0.069 

1256 

901 

1080 

1200 

1.083 

0.167 

6.002 

0.065 

1152 

849 

1020 

1130 

1.083 

0.223 

4.483 

0.022 

234 

213 

256 

284 

1.083 

0.223 

4.483 

0.044 

568 

426 

511 

568 

1.083 

0.223 

4.483 

0.066 

920 

629 

755 

839 

1.083 

0.223 

4.483 

0.091 

1200 

876 

1050 

1170 

1.083 

0.223 

4.483 

0.068 

900 

655 

786 

873 

1.083 

0.424 

2.360 

0.038 

140 

99 

119 

132 

1 .083 

0.424 

2.360 

0.075 

288 

198 

238 

264 

1.083 

0.424 

2.360 

0.099 

384 

261 

313 

348 

1.083 

0.424 

2.360 

0.099 

424 

260 

312 

346 

1.083 

0.424 

2.360 

0.091 

364 

239 

287 

319 

1.083 

0.424 

2.360 

0.088 

676 

231 

277 

308 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR RECTANGULAR MODEL 4 INCHES LONG 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in grams 

Cn.= '' *6 
m 

C =1 -8 

m 

C„,= 2.0 
m 

2.000 

0.173 

5.770 

0.031 

320 

323 

363 

404 

2.000 

0.173 

5.770 

0.038 

431 

389 

437 

486 

2.000 

0.173 

5.770 

0.045 

520 

468 

526 

584 

2.000 

0.173 

5.770 

0.045 

230 

236 

266 

295 

2.000 

0.173 

5.770 

0.014 

160 

149 

168 

186 

2.000 

0.345 

2.898 

0.037 

180 

141 

158 

176 

2.000 

0.345 

2.898 

0.055 

280 

210 

237 

263 

2.000 

0.345 

2.898 

0.074 

374 

283 

318 

354 

2.000 

0.345 

2.898 

0.075 

373 

287 

323 

359 

2.000 

0.345 

2.898 

0.068 

315 

259 

292 

324 

1 .500 

0.200 

4.992 

0.013 

no 

120 

135 

150 

1.500 

0.200 

4.992 

0.032 

245 

287 

323 

359 

1.500 

0.200 

4.992 

0.042 

375 

372 

419 

465 

1.500 

0.200 

4.992 

0.059 

554 

527 

593 

658 

1.500 

0.200 

4.992 

0.067 

716 

601 

676 

751 

1.500 

0.200 

4.992 

0.082 

800 

733 

824 

916 

1.500 

0.479 

2.086 

0.033 

64 

54 

60 

67 

1.500 

0.479 

2.086 

0.066 

132 

107 

121 

134 

1.500 

0.479 

2.086 

0.096 

200 

155 

174 

194 

1.500 

0.479 

2.086 

0.100 

210 

162 

182 

202 

1 .500 

0.479 

2.086 

0.094 

200 

152 

171 

190 

1.500 

0.479 

2.086 

0.087 

179 

141 

159 

176 

1 .083 

0.378 

2.642 

0.032 

112 

97 

109 

121 

1.083 

0.378 

2.642 

0.067 

230 

201 

226 

252 

1.083 

0.378 

2.642 

0.070 

230 

210 

236 

262 

1.083 

0.378 

2.642 

0.098 

320 

294 

330 

367 

1.083 

0.378 

2.642 

0.098 

330 

294 

330 

367 

1.083 

0.378 

2.642 

0.091 

300 

273 

307 

341 

1.083 

0.378 

2.642 

0.089 

310 

268 

302 

336 

1.083 

0.378 

2.642 

0.085 

304 

255 

287 

318 












J 








85 


COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR RECTANGULAR MODEL 4 INCHES LONG 
(Continued) 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in grams 

m 

00 

It 

E 

(_) 

C,„=2.0 

m 

1.083 

0.220 

4.543 

0.020 

140 

162 

182 

202 

1.083 

0.220 

4.543 

0.039 

300 

308 

346 

385 

1.083 

0.220 

4.543 

0.057 

460 

454 

510 

567 

1.083 

0.220 

4.543 

0.071 

600 

567 

638 

709 

1.083 

0.220 

4.543 

0.073 

625 

580 

652 

725 

1.083 

0.220 

4.543 

0.086 

750 

685 

771 

856 

1.083 

0.169 

5.910 

0.011 

124 

118 

133 

148 

1.083 

0.169 

5.910 

0.024 

265 

249 

280 

312 

1.083 

0.169 

5.910 

0.039 

457 

412 

463 

514 

1.083 

0.169 

5.910 

0.055 

694 

574 

646 

717 

1 .083 

0.169 

5.910 

0.069 

828 

720 

810 

900 

1 .083 

0.169 

6.910 

0.066 

806 

697 

784 

871 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR RECTANGULAR MODEL 2 INCHES LONG 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in grams 

m 

m 

m 

2.000 

0.180 

5.547 

0.006 

56 

49 

57 

65 

2.000 

0.180 

5.547 

0.014 

116 

108 

126 

144 

2.000 

0.180 

5.547 

0.015 

118 

112 

131 

150 

2.000 

0.180 

5.547 

0.025 

196 

184 

214 

245 

2.000 

0.180 

5.547 

0.035 

296 

262 

306 

350 

2.000 

0.180 

5.547 

0.043 

392 

320 

374 

427 

2.000 

0.180 

5.547 

0.060 

506 

447 

522 

596 

2.000 

0.180 

5.547 

0.061 

624 

460 

536 

613 

2.000 

0.344 

2.909 

0.021 

58 

61 

71 

81 

2.000 

0.344 

2.909 

0.043 

121 

125 

146 

.167 

2.000 

0.344 

2.909 

0.068 

183 

197 

230 

263 

2.000 

0.344 

2.909 

0.089 

244 

259 

302 

345 

2.000 

0.344 

2.909 

0.094 

296 

273 

318 

364 

2.000 

0.344 

2.909 

0.094 

266 

272 

317 

362 

2.000 

0.646 

1.547 

0.024 

5 

10 

12 

14 

2.000 

0.646 

1.547 

0.049 

23 

21 

24 

28 

2.000 

0.646 

1.547 

0.071 

35 

30 

35 

40 

2.000 

0.646 

1.547 

0.094 

48 

40 

47 

54 

2.000 

0.646 

1.547 

0.110 

52 

47 

54 

62 

2.000 

0.646 

1.547 

0.109 

58 

46 

54 

62 

1.500 

0.467 

2.140 

0.032 

43 

42 

50 

57 

1.500 

0.467 

2.140 

0.062 

94 

82 

95 

109 

1.500 

0.467 

2.140 

0.083 

132 

109 

127 

145 

1.500 

0.467 

2.140 

0.099 

142 

130 

152 

173 

1.500 

0.467 

2.140 

0.088 

110 

116 

135 

154 

1.500 

0.201 

4.975 

0.013 

78 

84 

98 

113 

1.500 

0.201 

4.975 

0.027 

182 

178 

207 

237 

1.500 

0.201 

4.975 

0.040 

290 

270 

315 

360 

1.500 

0.201 

4.975 

0.054 

462 

364 

425 

486 

1.500 

0.201 

4.975 

0.061 

568 

411 

480 

548 

1.500 

0.201 

4.975 

0.078 

640 

524 

612 

699 

1.500 

0.201 

4.975 

0.069 

524 

460 

536 

613 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR RECTANGULAR MODEL 2 INCHES LONG 
(Continued) 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 

Force in grams 

m 

m 

m 

1.083 

0.171 

5.862 

0.011 

82 

86 

101 

115 

1.083 

0.171 

5.862 

0.022 

184 

174 

203 

231 

1.083 

0.171 

5.862 

0.037 

364 

293 

342 

391 

0.083 

0.171 

5.862 

0.040 

392 

309 

361 

412 

1 .083 

0.171 

5.862 

0.057 

564 

446 

520 

594 

1.083 

0.171 

5.862 

0.070 

692 

549 

641 

732 

1.083 

0.171 

5.862 

0.068 

684 

529 

618 

706 

1.083 

0.171 

5.862 

0.064 

644 

500 

583 

667 

1.083 

0.221 

4.523 

0.019 

118 

112 

130 

149 

1.083 

0.221 

4.523 

0.043 

258 

257 

300 

343 

1.083 

0.221 

4.523 

0.067 

422 

397 

463 

529 

1.083 

0.221 

4.523 

0.093 

572 

556 

648 

741 

1.083 

0.221 

4.523 

0.071 

440 

422 

493 

563 

1.083 

0.377 

2.655 

0.031 

78 

70 

81 

93 

1.083 

0.377 

2.655 

0.059 

152 

135 

157 

179 

1.083 

0.377 

2.655 

0.093 

220 

211 

246 

281 

1.083 

0.377 

2.655 - 

0.094 

222 

214 

249 

285 

1.083 

0.377 

2.655 

0.077 

174 

175 

204 

233 













11 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR HALF CYLINDER MODEL (END FACING INCIDENT WAVE) 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in arams 

C =l-6 
m 

00 

1) 

E 

c_> 

Cn,= 2-0 

m 

2.000 

0.180 

5.552 

0.007 

82 

83 

93 

104 

2.000 

0.180 

5.552 

0.016 

206 

179 

202 

224 

2.000 

0.180 

5.552 

0.024 

344 

277 

311 

346 

2.000 

0.180 

5.552 

0.034 

510 

386 

434 

482 

2.000 

0.180 

5.552 

0.047 

720 

545 

613 

682 

2.000 

0.180 

5.552 

0.055 

800 

638 

718 

798 

2.000 

0.180 

5.552 

0.069 

988 

790 

889 

987 

2.000 

0.351 

2.850 

0.020 

88 

84 

94 

105 

2.000 

0.351 

2.850 

0.042 

184 

179 

202 

224 

2.000 

0.351 

2.850 

0.064 

286 

270 

304 

337 

2.000 

0.351 

2.850 

0.084 

398 

355 

399 

443 

2.000 

0.351 

2.850 

0.097 

516 

411 

462 

514 

2.000 

0.351 

2.850 

0.091 

448 

383 

431 

478 

2.000 

0.653 

1 .532 

0.021 

20 

13 

14 

16 

2.000 

0.653 

1.532 

0.047 

36 

29 

32 

36 

2.000 

0.653 

1.532 

0.070 

56 

43 

48 

54 

2.000 

0.653 

1 .532 

0.093 

76 

57 

64 

71 

2.000 

0.653 

1.532 

0.180 

86 

66 

74 

82 

2.000 

0.653 

1.532 

0.100 

86 

61 

68 

76 

1.500 

0.474 

2.108 

0.033 

72 

62 

69 

77 

1.500 

0.474 

2.108 

0.067 

141 

125 

141 

156 

1.500 

0.474 

2.108 

0.092 

212 

173 

194 

216 

1.500 

0.474 

2.108 

0.098 

220 

183 

206 

229 

1.500 

0.474 

2.108 

0.087 

162 

162 

183 

203 

1.500 

0.200 

5.006 

0.011 

127 

113 

127 

141 

1.500 

0.200 

5.006 

0.025 

288 

254 

286 

318 

1.500 

0.200 

5.006 

0.040 

480 

416 

468 

520 

1.500 

0.200 

5.006 

0.055 

724 

566 

637 

708 

1.500 

0.200 

5.006 

0.070 

984 

717 

807 

897 

1.500 

0.200 

5.006 

0.080 

1072 

827 

931 

1030 

1.500 

0.200 

5.006 

0.068 

864 

698 

785 

873 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR HALF CYLINDER MODEL (END FACING INCIDENT WAVE) 

(Continued) 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in grams 

m 

C =1.8 
m 

C„=2.0 

m 

1.083 

0.169 

5.926 

0.012 

128 

150 

169 

188 

1.083 

0.169 

5.926 

0.024 

320 

285 

321 

357 

1 .083 

0.169 

5.926 

0.038 

560 

464 

522 

580 

1.083 

0.169 

5.926 

0.059 

904 

714 

803 

892 

1.083 

0.169 

5.926 

0.069 

1152 

826 

930 

1030 

1.083 

0.169 

5.926 

0.062 

1096 

751 

845 

939 

1.083 

0.222 

4.512 

0.018 

172 

162 

183 

203 

1.083 

0.222 

4.512 

0.041 

406 

369 

415 

461 

1.083 

0.222 

4.512 

0.064 

684 

581 

653 

726 

1.083 

0.222 

4.512 

0.088 

952 

797 

896 

996 

1.083 

0.222 

4.512 

0.067 

720 

601 

676 

751 

1.083 

0.382 

2.619 

0.032 

104 

104 

117 

130 

1.083 

0.382 

2.619 

0.065 

220 

212 

239 

266 

1.083 

0.382 

2.619 

0.096 

336 

314 

353 

393 

1.083 

0.382 

2.619 

0.091 

300 

297 

334 

371 

1.083 

0.382 

2.619 

0.086 

288 

282 

317 

352 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR HALF CYLINDER MODEL (SIDE FACING INCIDENT WAVE) 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in grams 

C =1.6 
m 

O 

3 

II 

00 

o 

C\J 

II 

E 

o 

2.000 

0.179 

5.578 

0.008 

80 

89 

100 

111 

2.000 

0.179 

5.578 

0.016 

204 

179 

202 

224 

2.000 

0.179 

5.578 

0.025 

340 

284 

320 

356 

2.000 

0.179 

5.578 

0.034 

460 

390 

438 

487 

2.000 

0.179 

5.578 

0.045 

600 

515 

580 

644 

2.000 

0.179 

5.578 

0.056 

728 

643 

723 

804 

2.000 

0.179 

5.578 

0.066 

892 

754 

849 

943 

2.000 

0.344 

2.907 

0.021 

78 

91 

102 

113 

2.000 

0.344 

2.907 

0.043 

170 

189 

212 

236 

2.000 

0.344 

2.907 

0.065 

272 

286 

322 

358 

2.000 

0.344 

2.907 

0.087 

380 

382 

430 

477 

2.000 

0.344 

2.907 

0.086 

400 

377 

425 

472 

2.000 

0.344 

2.907 

0.096 

484 

423 

475 

528 

2.000 

0.344 

2.907 

0.095 

448 

417 

469 

521 

2.000 

0.653 

1 .532 

0.021 

18 

13 

14 

16 

2.000 

0.653 

1 .532 

0.021 

34 

27 

30 

34 

2.000 

0.653 

1.532 

0.021 

49 

38 

43 

47 

2.000 

0.653 

1 .532 

0.021 

69 

50 

56 

62 

2.000 

0.653 

1.532 

0.021 

74 

58 

65 

72 

2.000 

0.653 

1 .532 

0.021 

68 

55 

62 

69 

1.500 

0.473 

2.113 

0.030 

61 

57 

64 

71 

1.500 

0.473 

2.113 

0.030 

123 

118 

133 

148 

1.500 

0.473 

2.113 

0.030 

224 

173 

195 

217 

1.500 

0.473 

2.113 

0.030 

238 

169 

190 

212 

1.500 

0.473 

2.113 

0.030 

234 

173 

195 

217 

1 .500 

0.473 

2.113 

0.030 

170 

168 

189 

210 

1.500 

0.192 

5.206 

0.010 

121 

no 

123 

137 

1.500 

0.192 

5.206 

0.024 

318 

255 

287 

319 

1.500 

0.192 

5.206 

0.037 

492 

400 

450 

500 

1.500 

0.192 

5.206 

0.050 

720 

537 

604 

671 

1.500 

0.192 

5.206 

0.067 

900 

721 

811 

901 

1.500 

0.192 

5.206 

0.078 

1000 

830 

934 

1004 
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COMPARISON OF MEASURED AND COMPUTED HORIZONTAL 
FORCES FOR HALF CYLINDER MODEL (SIDE FACING INCIDENT WAVE) 

(Continued) 


d in 
feet 

d/L 

L/d 

H/L 

Measured 
Horizontal 
Force in 
grams 

Computed Horizontal 
Force in grams 

C =1.6 
m 

CO 

11 

E 

o 

C =2.0 
m 

1.083 

0.164 

6.088 

0.010 

151 

123 

138 

153 

1.083 

0.164 

6.088 

0.021 

370 

259 

291 

323 

1.083 

0.164 

6.088 

0.037 

604 

450 

506 

563 

1.083 

0.164 

6.088 

0.051 

860 

623 

701 

779 

1.083 

0.164 

6.088 

0.066 

964 

811 

912 

1010 

1.083 

0.164 

6.088 

0.061 

1008 

748 

841 

935 

1 .083 

0.164 

6.088 

0.057 

944 

703 

791 

879 

1.083 

0.207 

4.836 

0.017 

174 

164 

184 

205 

1.083 

0.207 

4.836 

0.040 

428 

393 

442 

491 

1.083 

0.207 

4.836 

0.061 

672 

600 

675 

750 

1 .083 

0.207 

4.836 

0.080 

916 

779 

876 

973 

1.083 

0.207 

4.836 

0.069 

700 

671 

754 

838 

1.083 

0.367 

2.725 

0.027 

100 

97 

110 

122 

1.083 

0.367 

2.725 

0.063 

260 

226 

255 

283 

1.083 

0.367 

2.725 

0.093 

406 

333 

375 

417 

1.083 

0.367 

2.725 

0.085 

402 

307 

345 

383 
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APPENDIX V 

SAMPLE COMPUTATION OF HORIZONTAL FORCE 

The following is a sample computation for a submerged tank 
utilizing the dimensionless force plots presented in Appendix II. 

A check on the answer obtained by the use of these curves is made 
by computing the force utilizing the method set forth by Reid 
and Bretschneider (15). 

Given: 

Wave; Period = T = 10 seconds 

Height = H = 10 feet 

Water depth = d = 150 feet 

Model: Height = H =40 feet 

m 

Width = W = 60 feet 
m 

Length = L =70 feet 
m 

Lq = 5.12 T^.(23) 

= 512 feet 

d/L^ = 150/512 = 0.293 

d/L = 0.3058 [from Appendix I of Wiegel (16)] 

L = 491 

L/d = 3.27 

H/L = 10/491 = 0.02 

''dim 


= 14.3 


[from Appendix II, Fig. 20 ] 



I 


1 


I 


•ii' 
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max 


DIM 


14.3 


a3 

Y A 

H 

d 

2 

*^max 

a3 

Y A 

H 

d 

2 

H 2 
m 


L 

m 



70 


11 .850 


•^max - —- 

= (14.3)(11850)(62.4)(10) 

(2T(150) 

= 352,000 LB 

Check; 

An inertial coefficient of 1.8 will be used in Eq. 40. 

ttL, 


(Fu)„,v = K Y H sin (-j^) 

H max m m m L 


K = 0.30 (from Fig. 41 or from Eq. 37] 

(FH)max = (1.8) (40) (60) (0.30) (62.4) (10) (sin(^^||^)) 

= (1.8)(40)(60)(0.30)(62.4)(10)(sin(0.448)) 
= 347,000 LBS 


(40) 
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ILO 


iT) 


_L_ 

lO 

o 


d/T-^Cf t./sec‘^) 

FI6.4I-PRESSURE FACTOR [AFTER REID AND 
BRETSCHNEIDER (15)] 
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APPENDIX VI 

COMPUTER PROGRAM FOR WAVE PARAMETERS AND INERTIAL WAVE FORCES 
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LINES = PARAMETER TO CONTROL CARRIAGE FOR PAGE CHANGE 
HORZFW = measured HORIZONTAL FORCE (IN GRAMS) IN DIRECTION OF WAVES 

HORZFO = MEASURED HORIZONTAL FORCE (IN GRAMS) IN OPPOSITE DIRECT FROM WAVE 

VERTFD = measured VERTICAL FORCE (IN GRAMS) DOWN 

VERTFU = MEASURED VERTICAL FORCE (IN GRAMS) UP 







HORZW = DIMENSIONLESS HORIZONTAL FORCE IN DIRECTION OF WAVES 

HORZO = DIMENSIONLESS HORIZONTAL FORCE IN OPPOSITE DIRECTION FROM WAVES 

VERTD = DIMENSIONLESS VERTICAL FORCE DOWN 

VERTU = DIMENSIONLESS VERTICAL FORCE UP 

GAMMA = unit WEIGHT OF WATER 
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FUNCTION SINH(X) 730 
SINH=(EXP(X)-(l./EXP(X)))/2. 740 
return 750 
END 760 
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APPEi'iniX VII 


COMPUTER PROGRAM FOR OBJECTS OF IRREGULAR CROSS SECTIOM 
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THE NEXT FOUR CAROS ARE FOR ALL MODELS EXCEPT FLAT PLATE 
THEY CONVERT FORCE IN GRAMS TO DIMENSIONLESS FORCE. 
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FUNCTION SINH(X) 73O 
SINH=(EXP(X)-ll./EXP(X)))/2. 74 C 
RETURN 750 
END 760 
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